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Abstract. 

Completely integrable Hamiltonian systems look rather promising for controllability since their 
first integrals are stable under an internal evolution, and one may hope to find a perturbation of a 
Hamiltonian which drives the first integrals at will. Action-angle coordinates are most convenient for 
this purpose. Written with respect to these coordinates, a Hamiltonian and first integrals of a (time- 
dependent) completely integrable system depend only on action variables. We introduce a suitable 
perturbation of an internal Hamiltonian by a term containing time-dependent parameters (control 
fields) so that an evolution of action variables is nothing else than a holonomy displacement along 
a curve in a parameter space. Therefore, one can determine this evolution in full by an appropriate 
choice of control fields. Similar holonomy controllability of finite level quantum systems is of special 
interest in connection with quantum computation. We provide geometric quantization of a time- 
dependent completely integrable Hamiltonian system in action-angle variables. Its Hamiltonian 
and first integrals have time-independent countable spectra. The holonomy control operator in this 
countable level quantum system is constructed. 

1 Introduction 

A time-dependent Hamiltonian system of m degrees of freedom is called a completely inte- 
grable system (henceforth CIS) if there exist m independent first integrals in involution. We 
show that such a system admits the action-angle coordinates around any regular instantly 
compact invariant manifold. Written relative to these coordinates, its Hamiltonian and first 
integrals are functions only of action coordinates. In comparison with perturbations in the 
KAM theorem, we consider perturbations of a Hamiltonian of a CIS by the term which 
contains time-dependent parameters. A generic Hamiltonian of a mechanical system with 
time-dependent parameters includes a term which is linear in momenta and the temporal 
derivatives of parameter functions [15, 11, 2^| . Then the corresponding evolution operator 



of action variables depends only on a trajectory of the parameter functions in a parameter 
space. Therefore, it is determined completely by an appropriate choice of these parameter 
functions (control fields), and plays the role of a holonomy control operator. 

At present, holonomy control operators in quantum systems attract special attention 
in connection with quantum computation, based on the generalization of Berry's phase by 
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means of driving a finite level degenerate eigenstate of a Hamiltonian over the parameter 
manifold || ^U], [28|]. Information is encoded in this degenerate state. Bearing in mind 
this application, we quantize a time-dependent CIS in question, and construct the quantum 
holonomy control operator in this system. 

An essential simplification is that, given action- angle variables, a time-dependent CIS 
can be quantized as an autonomous one because its Hamiltonian and first integrals relative 
to these coordinates are time- independent. Of course, the choice of action-angle coordinates 
by no means is unique. They are subject to canonical transformations. Therefore, we employ 



the geometric quantization technique || ^5|, which remains equivalent under canonical 
transformations, but essentially depends on a choice of polarization |4], |2]J . 

Geometric quantization of an autonomous CIS has been studied with respect to polar- 
ization spanned by Hamiltonian vector fields of first integrals [|Tj|. The well-known Simms 
quantization of the harmonic oscillator is also of this type. The problem is that the associ- 
ated quantum algebra includes functions which are ill defined on the whole momentum phase 
space, and elements of the carrier space fail to be smooth sections of the quantum bundle. 
Indeed, written with respect to the action-angle variables, this quantum algebra consists of 
functions which are affine in angle coordinates. 

We choose a different polarization spanned by almost-Hamiltonian vector fields of angle 
variables. The associated quantum algebra A consists of smooth functions which are affine 
in action variables. This quantization is equivalent to geometric quantization of the cotan- 
gent bundle of the torus T m with respect to the vertical polarization. This polarization is 
known to lead to Schrodinger quantization. We show that A possesses a set of inequivalent 
representations in the separable pre-Hilbert space C°°(T m ) of smooth complex functions on 
T m . In particular, the action operators read 

I k = -id k - A fc , (1) 

where are real numbers which specify different representations of A. By virtue of the 
multidimensional Fourier theorem, an orthonormal basis for C°°(T m ) consists of functions 

VV) = exp[i(n r r )], (n r ) = (m, . . . , n m ) G Z m , (2) 

where <f) 1 are cyclic coordinates on T m . With respect to this basis, the action operators (|l|) 
are written as countable diagonal matrices 

4VV) = ( n k ~ ^k)lp(n r )- (3) 
Given the representation (|1|), any polynomial Hamiltonian H(Ik) of a CIS is uniquely quan- 
tized as a Hermitian element Tt{Ik) = W(ljt) of the enveloping algebra of A. It has the 
time-independent countable spectrum 

H{h)^(n r ) = £(„ r )VV)' E (n r ) = 7~L{n k - Afc), n k G (n r ). (4) 

Since Ik are diagonal, one can also quantize Hamiltonians H{Ij) which are analytic functions 
on R m . 

Quantum CISs look rather promising for quantum computation. Its Hamiltonian depends 
only on action variables and possesses a time-independent countable spectrum. Moreover, 
it can be made degenerate at will by appropriate canonical transformations. We construct 
the holonomy control operator in this countable level quantum system. Note that one used 
to study controllability of finite level quantum systems |H, 3, ^4}|. 
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2 Classical time-dependent completely integrable systems 



In order to introduce action-angle variables for a time-dependent CIS, we use the fact that 
a time- dependent Hamiltonian system of m degrees of freedom can be extended to an au- 
tonomous one of m + 1 degrees of freedom where the time is treated as a dynamic variable 
P, [7|, P^|. The classical theorem || [13j on action-angle coordinates around a regular compact 
invariant manifold can not be applied to this autonomous CIS since its invariant manifolds 
are never compact because of the time axis. Generalizing the above mentioned theorem, we 
first prove that there is a system of action-angle coordinates on an open neighbourhood U of 
a regular connected invariant manifold M of an autonomous CIS if Hamiltonian vector fields 
of first integrals on U are complete and the foliation of U by invariant manifolds is trivial. 
If M is a compact regular invariant manifold, these conditions always hold |TB[ . Afterwards, 



we show that, if a regular connected invariant manifold of a time-dependent CIS is compact 
at each instant, it is diffeomorphic to the product of the time axis R and an m-dimensional 
torus T m , and admits an open neighbourhood endowed with time-dependent action- angle 
coordinates (If, t, <ft), % — 1, . . . , m, where t is the Cartesian coordinate on R and <p l are cyclic 
coordinates on T m . 

Recall that the configuration space of a time- dependent mechanical system is a fibre 
bundle Q — > R over the time axis R equipped with the bundle coordinates (t,q k ), k = 
1, . . . ,m. The corresponding momentum phase space is the vertical cotangent bundle V*Q 
of Q — > R endowed with holonomic coordinates (t, q k ,Pk = q k ) 13 The cotangent bundle 
T*Q, coordinated by (q x ,P\) = (t,q k ,p ,p k ), plays a role of the homogeneous momentum 
phase space. It is provided with the canonical Liouville form H = p\dq x , the canonical 
symplectic form Q = dp\ A dq x , and the corresponding Poisson bracket 

{/, f'h = d x fd x f - d x fd x f, f, f e C°°(T*Q). (5) 

There is the one-dimensional trivial affine bundle 

C : T*Q -> V*Q. (6) 

Given its global section h, one can equip T*Q with the global fibre coordinate r = p — h. 
The fibre bundle (0) provides the vertical cotangent bundle V*Q with the canonical Poisson 
structure {, }y such that 

C*{/, f'}v = {Cf, Cf'h, V/, /' e C°°(V*Q), (7) 
{f,f}v = d k fd k f-d k fd k f. (8) 

A Hamiltonian of time- dependent mechanics is defined as a global section 

h:V*Q-*T*Q, p oh = -H(t,q j , Pj ), (9) 

of the affine bundle ( (((]) [IS, It yields the pull-back Hamiltonian form 

H = h*E = p k dq k - Hdt 

on V*Q. There exists a unique vector field 7# on V*Q such that 

lH \dt = l, lH \dH = 0, 

lH = d t + d k Hd k -d k Hd k . (10) 
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Its trajectories obey the Hamilton equation 

q k = d k n, p k = -d k n. (ii) 

A first integral of the Hamilton equation (pTTj) is defined as a smooth real function F on 
V*Q whose Lie derivative 

L lH F = lH \dF = d t F + {H,F} v 

along the vector field jh (|10D vanishes, i.e., F is constant on trajectories of 7#. A time- 
dependent Hamiltonian system (V*Q, H) is said to be completely integrable if the Hamilton 
equation ([II]) admits m first integrals F k which are in involution with respect to the Poisson 
bracket {, }y (§), and whose differentials dF k are linearly independent almost everywhere 
(i.e., the set of points where this condition fails is nowhere dense). One can associate to this 
CIS an autonomous CIS on T*Q as follows. 

Given a Hamiltonian h it is readily observed that 

H* = d t \(E-Ch*E)= Po + H (12) 

is a function on T*Q. Let us regard TC* as a Hamiltonian of an autonomous Hamiltonian 
system on the symplectic manifold (T*Q,Q) [[L7|]. Its Hamiltonian vector field 

7r = d t - d t Hd° + d k Hd k - d k Hd h (13) 

is projected onto the vector field 7^ ( |10D on V*Q so that 

C*(L 7H /) = {H\ Cfh, V/ e C°°(V*Q). 

An immediate consequence of this relation is the following. 

Proposition 1. (i) Given a time-dependent CIS (H;F k ) on V*Q, the Hamiltonian sys- 
tem (H*,C*Fk) 011 T*Q is a CIS. (ii) Let iV be a connected regular invariant manifold of 
(H;Fk). Then h(N) C T*Q is a connected regular invariant manifold of the autonomous 
CIS (H*, (*F k ). 

Hereafter, we assume that the vector field 7^ ( [T0| ) is complete. In this case, the Hamilton 
equation (|Tl]) admits a unique global solution through each point of the momentum phase 
space V*Q, and trajectories of -yn define a trivial bundle V*Q — > V t *Q over any fibre V*Q of 
V*Q —> R. Without loss of generality, we choose the fibre i : VqQ — > V*Q at t — 0. Since 
is an invariant manifold, the fibration 

£ : V*Q -> V*Q (14) 

also yields the fibration of iV onto N — N n V^Q such that = R x jV Q is a trivial bundle. 
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3 Time-dependent action-angle coordinates 



Let us introduce the action-angle coordinates around an invariant manifold N of a time- 
dependent CIS on V*Q by use of the action-angle coordinates around the invariant manifold 
h(N) of the autonomous CIS on T*Q in Proposition |l|. Since N and, consequently, h(N) are 
non-compact, we refer to the following. 

Proposition 2. Let M be a connected invariant manifold of an autonomous CIS {Fx}, 
A = 1, . . . , n, on a symplectic manifold (Z, Viz)- Let U be an open neighbourhood of M such 
that: (i) the differentials dF x are independent everywhere in U, (ii) the Hamiltonian vector 
fields $\ of the first integrals F\ on U are complete, and (iii) the submersion xF\ : U — > M n 
is a trivial bundle of invariant manifolds over a domain 7' C l n . Then U is isomorphic to 
the symplectic annulus 

W' = V'x (R n - m x T m ), (15) 
provided with the action-angle coordinates 

(/!,..., 4; x\...,x n - m -cf>\...,<t> m ) (16) 

such that the symplectic form on W reads 

Q z = dI a A dx a + dli A d(j)\ 

and the first integrals Fa depend only on the action coordinates la- 
Proof. In accordance with the well-known theorem ||, the invariant manifold M is diffeo- 
morphic to the product M n_m x T m , which is the group space of the quotient G = M n /Z m of 
the group lR n generated by Hamiltonian vector fields $a of first integrals Fx on M. Namely, 
M is provided with the group space coordinates (y x ) = (s a , tp l ) where ip % are linear functions 
of parameters s x along integral curves of the Hamiltonian vector fields $\ on U. Let (J\) 
be coordinates on V which are values of first integrals Fx- Let us choose a trivialization of 
the fibre bundle U — > V seen as a principal bundle with the structure group G. We fix its 
global section x- Since parameters s x are given up to a shift, let us provide each fibre Mj, 
J G V, with the group space coordinates (y x ) centred at the point x(J)- Then (J x ;y x ) are 
bundle coordinates on the annulus W (0) . The rest of the proof in Appendix A reduces to 
transformation of the coordinates ( Jx] y x ) to the desired coordinates 

L = Ja, Ii(Jj), X a = S a + S a (J X ), ct> l = ^ + S\J X ,S b ). (17) 



Of course, the action-angle coordinates (|TE|) by no means are unique. For instance, let 
T a be (n — m) arbitrary smooth functions on IR m . Let us consider the canonical coordinate 
transformation 

l'a = Ia-^a{Ij), 4=4, x' a = x\ ^=^ + 1^(1,). (18) 

Then (I' a , I' k ; x' a , <p' k ) are action-angle coordinates on the symplectic annulus which differs 
from W fll5| ) in a trivialization. 
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Let us apply Proposition § to the CISs in Proposition [T]. 



Proposition 3. Let N be a connected regular invariant manifold of a time- dependent CIS on 
V*Q, and let the image N of its projection £ ([14]) be compact. Then the invariant manifold 
fo(iV) of the associated autonomous CIS on T*Q has an open neighbourhood U obeying the 
condition of Proposition Q (see Appendix B). 

In accordance with Proposition ||], the open neighbourhood U of the invariant manifold 
h(N) of the autonomous CIS in Proposition |3] is isomorphic to the symplectic annulus 

W' = V x (R x T m ) (19) 

provided with the action-angle coordinates (Jo, . . . , I m ; t, 1 , . . . , m ) such that the symplectic 
form on W reads 

Q = dl A dt + dl k A 

By the construction in Proposition I = J = 7i* ( |T2"D and the corresponding angle 
coordinate is x° = t, while the first integrals Jt = (*Fk depend only on the action coordinates 

u 

Since the action coordinates Jj are independent of the coordinate Jo, the symplectic 
annulus W (|19[) inherits the fibration (I6T) which reads 



C:W'3 (J , h; t, <jj) -> (Ji, t, <?) e W = R x T m x V. (20) 

By the relation similar to (^), the product (|20| ) is provided with the Poisson structure 

{/, /V = VfdJ' - /, /' e C°°(M/). 

Therefore, one can regard IV as the momentum phase space of the time-dependent CIS in 
question around an invariant manifold N. 

It is readily observed that the Hamiltonian vector field of the autonomous Hamiltonian 
7i* = Io is 7t = d t , and so is its projection 7# (|I0D on VF. Consequently, the Hamilton 
equation ( pTT]) of a time-dependent CIS with respect to the action-angle coordinates take the 
form Jj = 0, cff = 0. Hence, (If, t, (p 1 ) are the initial data coordinates. One can introduce such 
coordinates as follows. Given the fibration £ flUD , let us provide N xV C VqQ in Proposition 
H with the action-angle coordinates (If,~(j)) for the CIS {io-^fc} on the symplectic leaf VqQ. 
Then, it is readily observed that (If t, 4> ) are time-dependent action-angle coordinates on W 



(|20|) such that the Hamiltonian H.(Ij) of a time-dependent CIS relative to these coordinates 
vanishes, i.e., 7i* = Iq. Using the canonical transformations (|IS|), one can obtain different 
time-dependent action- angle coordinates. In particular, given a smooth function 7i on M. m , 
one can endow W with the action-angle coordinates 

I =7 -H(I J ), h=li, <p l = ^ +td l H(I j ) 
such that TC(I{) is a Hamiltonian of time-dependent CIS on W. 
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4 The classical control operator 

A generic momentum phase space of a Hamiltonian system with time-dependent parameters 
is a composite fibre bundle II — > E — >• R, where II — > S is a symplectic bundle and X — > R 



is a parameter bundle whose sections are parameter functions [11, Pj|, E3L E7l. Here, we 



assume that all bundles are trivial and, moreover, their trivializations hold fixed. Then the 
momentum phase space of a Hamiltonian system with time-dependent parameters on the 
Poisson manifold W ( p0| ) is the product 

U = SxW = Vx(RxSxT m )^RxS^R, (21) 

equipped with the coordinates (/&; t, <7 a , fc ). It is convenient to suppose for a time that 
parameters are also dynamic variables. The momentum phase space of such a system is n' = 
T*S x W, coordinated by (Ik]t,a a ,p a ,(j) k ). The dynamics of a time-dependent mechanical 
system on the momentum phase space n' is characterized by a Hamiltonian form 

#s = Pada a + I k d(j) k - Hv(t, a^pp, J„ <P)dt 

fts = Pa T? + I k (A k t + T?A*) + H, (22) 

where T = (Tf) is a connection on the parameter bundle R x S —>■ R and A = (A k , A*) is a 
connection onlxSxT m ^lx5 [[TT], [23[ . 



Bearing in mind that a a are parameters, one should choose the Hamiltonian TC^ (22) to 
be affine in momenta p a . Furthermore, in order to describe a Hamiltonian system with a 
fixed parameter function <j a = £, a {t), one defines the connection V such that 

v r e = o, r?(t,^(t)) = d t c. 

Then the pull-back 

H 6 = = W - (I*[A*(f , e, V) + Aj(f , ^)d t C] + H(t, Ij, P))dt 
is a Hamiltonian form on the Poisson manifold W (^). Let us put 

H = H — I k A t , 

where TC(Ii) is a Hamiltonian of the original CIS on W (|20|). Then the Hamiltonian form 

ff € = - Hzdt = I k d<p k - [J fc A*(t, tf, <tP)d t C + W&OJdi (23) 

describes a perturbed time-dependent CIS on the Poisson manifold W ( pfj|) . The correspond- 
ing Hamilton equation reads 

a t / fc = -dkAiijdtt", d t <p k = d k n + A k a d t c. (24) 

In order to make the term 

A = I k A%C (25) 

in the perturbed Hamiltonian (|23| ) a control operator, let us assume that the coefficients 
A^ of the connection A are independent of time. Then, in view of the trivialization (21), 
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its part (A*) can be seen as a connection on the fibre bundle S x T m — > S. Let us choose 
the initial data action-angle variables. Then the internal Hamiltonian TC of the original CIS 



vanishes, and the Hamilton equation (|24 ) takes the form 



(26) 
(27) 



This is the control equation as follows. 

Any smooth complex function on the product RxT m is represented by a multidimensional 
Fourier series of functions ipi nr ) @ on T m with coefficients the smooth functions on R. Let 
us rewrite the equation ( |27D as the countable system of equations 

d t ip(n r ) = iip(n r )nid t (j)' = iVV0 n * A a(£^VVr))dt<r, rii G (n r ), 

for functions VW) (0). Let 

K^J2K(m r )(^)Amr) (28) 
(m r ) 

be the Fourier series for A l a . Since T/VO'^m,.) = ^{n r +m r )i we obtain a countable system of 
linear ordinary differential equations 

= E[^Si(^)^r]^w, (29) 



M 



Or) 



a(n r ) 



(mr-)' 



(30) 



(m r +k r =n r ) 



with time-dependent coefficients [iM^f\{^)dti a \. Its solution with the initial data ip l (0) 
can be written as the formal time-ordered exponential 



U(t) = Texp 



i J M a (e(t'))d t Cdt' 
o 



Texp 



/ / M a (a^)da c 
«[o,t]) 



(31) 



where M Q denotes the matrix with time-dependent entries (|30|) [12], [1{J. A glance at the evo- 
lution operator U(t) ([31]) shows that solutions of the equations ( p9| ) are functions ^(n r ) (£(£)) 
of a point £(t) of the curve £ : R — > S 1 in the parameter space S. 

Substituting this solution into the equation (|26|), we obtain the system of I ordinary linear 
differential equations with time-dependent coefficients: 

d t I t = ~[L^(t),i; {nr) m))9tC}I k , L k m = ftAj 



Its solution is given by the time-ordered exponential 

hit) = U(t)tl h (0), 
t 

'L a (e(t)^ {nr) (at)))dtCdt' 



U(t) = Texp 



Texp 



L^a 13 ,ip{n r) (a))d. 



a 



f([o,t]) 
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where L a denotes the matrix with time-dependent entities L k ai . This solution is a functions 
of a point £ if) of the curve £ : R — > S in the parameter space S. 

It follows that, if the holonomy group of the connection (A^) on the fibre bundle S x T m — > 
S is the whole group GL(m, M), one can obtain any desired trajectory of action variables Ii 
and, consequently, of first integrals by an appropriate choice of control functions 



5 Quantum completely integrable systems 

In order to quantize a time-dependent CIS on the Poisson manifold (W, {, }vy), one may follow 
the general procedure of instantwise geometric quantization of time-dependent Hamiltonian 
systems in [[lCj. As was mentioned above, it can however be quantized as an autonomous 



CIS on the symplectic annulus 

P = V xT r 

equipped with fixed action-angle coordinates and provided with the symplectic form 

tt P = dh A d<f>\ (32) 



In accordance with the standard geometric quantization procedure [^5], ^6| , because the 
symplectic form Qp (^) is exact, the prequantum bundle is defined as a trivial complex line 
bundle C over P. Since the action-angle coordinates are canonical for the symplectic form 
fl3"2]), the prequantum bundle C need no metaplectic correction. Let its trivialization 

C = P x C (33) 

hold fixed. Any other trivialization leads to equivalent quantization of P. Given the asso- 
ciated bundle coordinates 4> k , c), c G C, on C (|33"D, one can treat its sections as smooth 
complex functions on P. 

The Konstant-Souriau prequantization formula associates to each smooth real function 
/ G C°°(P) on P the first order differential operator 

f = -iVt f +f (34) 

on sections of C, where *&f = d k fdk — dkfd h is the Hamiltonian vector field of / and V is 
the covariant differential with respect to a suitable [/(l)-principal connection on C. This 
connection preserves the Hermitian metric g(c, c') = cc' on C, and its curvature form obeys 
the prequantization condition R = iflp. This connection reads 

A = A + icl k d(f) k <g> d c , (35) 

where A is a flat [/(l)-principal connection on C — > P. The equivalence classes of flat 
principal connections on C are indexed by the set M m /Z m of homomorphisms of the de 
Rham cohomology group 

H\P) = H\T m ) = R m 
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of P to the cycle group U(l) f|. We choose their representatives of the form 

A [(\ k )\ = dl k ®d k + d(j) k ® (d k + i\ k cd c ), X k e [0, 1). 
Then the connection up to gauge transformations reads 

A[(X k )] = dl k ®d k + d<p k <g> (d k + i{I k + \ k )cd c ). (36) 



For the sake of simplicity, we will assume that the numbers X k in the expression (p6|) belong 
to K, but will bear in mind that connections A[(A/-)] and A[(Aj.)] with X k — \' k G Z are gauge 
conjugated. Given a connection (0), the prequantization operators (|34]) read 



f = -i$ f + (f-(I k + \ k )d k f). (37) 

Let us choose the above mentioned angle polarization W which is the vertical tangent 
bundle of the fibration it : P — > T m , and is spanned by the vectors It is readily observed 
that the corresponding quantum algebra A C C°°(P) consists of affine functions 

f = a k (^)I k + b(^) (38) 

of action coordinates I k . The carrier space of its representation by operators ([J7|) is defined 
as the space E of sections p of the prequantum bundle C of compact support which obey 
the condition V^p = for any Hamiltonian vector field d subordinate to the distribution 
Vn. This condition reads 

d k fd k p = 0, V/GC°°(T m ). 

It follows that elements of E are independent of action variables and, consequently, fail to 
be of compact support, unless p = 0. This well-known problem of Schrodinger geometric 
quantization is solved as follows || [Tofl ■ 

Let us fix a slice '■ T m T m x V. Let CV = V^C be the pull-back of the prequantum 
bundle C (^) over the torus T m . It is a trivial complex line bundle Ct = T m x C provided 
with the pull-back Hermitian metric g(c, c') = cc'. Its sections are smooth complex functions 
on T m . Let 

A = i* T A = d(j) k ® (d k + i(I k + \ k )cd c ) 

be the pull-back of the connection A onto Ct- Let Dbea metalinear bundle of complex 
half-forms on the torus T m . It admits the canonical lift of any vector field r on T m , and the 
corresponding Lie derivative of its sections reads 

L r = r k d k + l -d k r k . 

Let us consider the tensor product 

Y = C T ®V ^T m . (39) 

Since the Hamiltonian vector fields 

y = a k d k - (I r d k a r + d k b)d k 
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of functions / fl3§D are projectable onto T m , one can associate to each element / of the 
quantum algebra A the first order differential operator 



/ = HV^. + /) <8> Id + Id ® L^ f = -ia k d k - l -d k a k - a k \ k + b (40) 



on sections of Y. A direct computation shows that the operators ( |4"0"[ ) obey the Dirac 
condition 

[/, /'] = -i{U'}- 

Sections pr of the quantum bundle Y — > T m constitute a pre-Hilbert space Et with 
respect to the non-degenerate Hermitian form 



(Pt\p' t ) = (^A J PtPt, Pt,Pt^E t . 



Then it is readily observed that / (|40D are Hermitian operators in Et- In particular, the 
action operators take the form ([1]). 

Of course, the above quantization depends on the choice of a connection A[(Afc)] (^6|) 
and a metalinear bundle T>. The latter need not be trivial. If T> is trivial, sections of the 
quantum bundle Y — > T m (^) obey the transformation rule 

p T {4> k + 2tt) = p T ^ k ) 

for all indices k. They are naturally complex smooth functions on T m . In this case, Et is 
the above mentioned pre-Hilbert space C°°(T m ) of complex smooth functions on T m whose 
basis consists of functions (0). The action operators I ([I]) with respect to this basis are 
represented by countable diagonal matrices (0), while functions a(<f>) are decomposed into 
the pull-back functions ip( nr ) which act on C°°(T m ) by multiplications 

VV)W r ) = VV)^«) = i>{n r +n' r )- (41) 

If D is a non-trivial metalinear bundle, sections of the quantum bundle Y — * T m (|39|) 
obey the transformation rule 

p T (<P j +2n) = -pr(p) (42) 

for some indices j. In this case, the orthonormal basis of the pre-Hilbert space Et can be 
represented by double-valued complex functions 

= ex P[*( n i0 l + ( n j + ^)<P J )} ( 43 ) 

on T m . They are eigenvectors 

of the operators I k ([]]), and the functions a(0) act on the basis ( fi5j) by the above law 
(|4l|). It follows that the representation of A determined by the connection A[(Afe)] (^6|) 
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in the space of sections ( f42|) of a non-trivial quantum bundle Y (|39|) is equivalent to its 
representation determined by the connection A[(A», \j — |)] in the space C°°(T m ) of smooth 
complex functions on T m . 

Therefore, one can restrict the study of representations of the quantum algebra A to its 
representations in C°°(T m ) associated to different connections (j36|). These representations 
are inequivalent, unless \ k — A' fc G Z for all indices k. 

Now, in order to quantize the Poisson manifold (W, {, }w), one can simply replace func- 
tions on T m with those on I x T m [pi], p5fl . Let us choose the angle polarization of W 
spanned by the vectors d k . The corresponding quantum algebra Aw C C°°{W) consists of 
affine functions 

f = a k (t^)I k + b(t^) (44) 

of action coordinates I k , represented by the operators (|4"0D in the space C°°(R x T m ) of 
smooth complex functions on 1 x T m . This space is provided with the structure of the 
pre-Hilbert C°°(R)-module with respect to the non-degenerate C°°(R)-bilinear form 

/ 1 \ m „ 

= {^) J W , 1>, # G x T m ). 

Its basis consists of the pull-backs onto R x T m of the functions ip( nr ) ©■ 

Since the Poisson structure (|2l|) defines no dynamics on the momentum phase space 
W (0), we should quantize the homogeneous momentum phase space W (|T9| ) in order 
to describe evolution of a quantum time-dependent CIS. Following the general scheme in 
jlOl , one can provide the relevant geometric quantization of the symplectic annulus 
(W',fl'). The corresponding quantum algebra Aw' G C°°(W') consists of affine functions 

f = a\t^)h + b(t^) 

of action coordinates I\. It suffices to consider its subalgebra consisting of the elements / 
and Iq + f for all / G Aw (@)- They are represented by the operators / (fiDD and To = —id t 
in the pre-Hilbert module C°°(R x T m ). If a Hamiltonian Ti{Ij) of the time-dependent CIS is 
a polynomial (or analytic) function in action variables, the Hamiltonian 7i* of the associated 
autonomous CIS is quantized as 

H* = -idt + Hifj). 
Then we obtain the Schrodinger equation 

H*i) = -idttp + H(-id k - X k )tp = 0, tp G C°°(R x T m ). 
Its solutions are the series 

i> = B M exp[-i-B(n r )#(rv)) B {n r ) G C, 

(n r ) 

where Er^ are the eigenvalues @) of the Hamiltonian Tt. 
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6 The quantum control operator 



In comparison with the classical control operator in Section 4, we will construct the quantum 
control operator which preserves the eigenvalues of a Hamiltonian, and acts in its degenerate 
eigenspaces. For instance, let us choose action-angle coordinates such that a Hamiltonian 
7i of a CIS equals Ji, and it is independent of other action variables I a (a, b, c — 2, . . . , m). 
It is quantized as Ti = l\. Its eigenvalues are countably degenerate. Let us consider the 
perturbed Hamiltonian (|23| ) where the perturbation term A ( p5| ) depends only on the 
action-angle coordinates with the indices a, b, c = 2, . . . , m, i.e., 



The perturbation term 



A = A*(e\0 6 )^/ a 



of this Hamiltonian is an element of the quantum algebra Aw, and is quantized by the 
operator 



£ i(k a + -m a -\ a )A« {ma) (e)], 

(m a +k a =n a ) 



j3(ni,n a ) 



where the Fourier series decomposition (|28| ) is used. 

Since the operators A and Ti mutually commute, the corresponding quantum evolution 
operator reduces to the product 



Texp 



Hedt' 









t 




t 


U 1 {t)oU 2 {t)=Texp 


-i J Hdt' 


o Texp 


-i J Adt' 













(45) 



The first factor in this product is the dynamic evolution operator of the quantum CIS. It 
reads 

Ul(t)lp( ni , na ) = exp[-z(m - Ai)*]V>(ni,ne)- ( 46 ) 

Its eigenvalues are countably degenerate. Recall that the operator (|46| ) acts in the pre- 
Hilbert module C°°(R x T m ). Its eigenvalues are smooth complex functions on R, and its 
eigenspaces are C°°(R)-submodules of C°°(R x T m ) of countable rank. 
The second factor in the product (flop is 



U 2 (t) = Texp 



-i j Ap(e(t'))d t fdt' 
o 



Texp 



-/ / A p (a a )da p 

am) 



(47) 



It acts as a matrix of countable rank in the eigenspaces of the internal Hamiltonian Ti. Its 
eigenspace corresponding to the eigenvalue {rij — Xj) is the pre-Hilbert C°°(R)-submodule of 
C°°(R x T m ) whose orthonormal basis is made up by functions ip( ni ,n a ) for all collections of 
integers {n a ). 
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A glance at the expression fl47|) shows that, in fact, the operator U2(t) depends on the 
curve £([0,1]) C S in the parameter space S. One can treat it as an operator of parallel 
displacement with respect to a connection in the C°°(E)-module of smooth complex functions 



on E x T m along the curve £ |TT|, [Ll| [ITJ. For instance, if £([0, 1]) is a loop in S, the operator 



U2 (0) is the geometric Berry factor. In this case, one can think of U2 as being a holonomy 
control operator. 

It should be emphasized that, in comparison with operators usually studied \% [| [| [20] 



24 , 28fl , the operator (|47D acts in a countable level quantum system. Of course, the problem 



arises if such a system admits complete controllability and if the holonomy control operator 
(|47| ) can provide this controllability. This problem will be studied elsewhere. 



7 Appendix A 

Let us complete the proof of |2|. Since Mj are Lagrangian manifolds, the symplectic form Vtz 
on W is given relative to the bundle coordinates ( J x ; y x ) by the expression 

Vl z = tt Q(3 dJ a A dJp + QpdJ a A dyP. (48) 

By the very definition of coordinates {y x ), the Hamiltonian vector fields $\ of first integrals 
take the coordinate form $ x = $"( J^)d a . Moreover, since the cyclic group S 1 can not act 
transitively on M, we have 

0a = da + KWdi, $i=^(Jx)d k . (49) 

The Hamiltonian vector fields i9\ obey the relations 

$ x \n z = -dj x , npl = 5 a x . (50) 

It follows that VL'p is a non-degenerate matrix and = (fi -1 )", i.e., the functions depend 
only on coordinates J\. A substitution of ( f4"9"D into (^) results in the equalities 

n a b = s a b , tfjni = 0, (51) 
= 51 tffng = 0. (52) 

The first of the equalities d5"2|) shows that the matrix fi^. is non-degenerate, and so is the 
matrix d^. Then the second one gives Q% = 0. 

By virtue of the well-known Kiinneth formula for the de Rham cohomology of a product 
of manifolds, the closed form Qz (fH) ori W ( pT5|) is exact, i.e., fi^ = dE where 5 reads 



5 = ~ a (J x , y x )dJ a + EiiJjdtp* + <9 Q $(J A , y x )dy a , 

where $ is a function on W. Taken up to an exact form, 5 is brought into the form 

3 = ~' a ( J x , y x )dJ a + 3,(J A ) dtp*. (53) 

Owing to the fact that components of cG = Qz are independent of y x and obey the equalities 
dH) _ (HD) we obtain the following. 
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(i) fi? = — (9jS /a + <9 a Hj = 0. It follows that d0 a is independent of p\ i.e., E' a is affine in 
tp l and, consequently, is independent of (p l since are cyclic coordinate. Hence, <9 a Sj = 0, 
i.e., Sj is a function only of coordinates Jj. 

(ii) f2f = —diE' k + <9 fc Sj. Similarly to item (i), one shows that E' k is independent of ip 1 
and Q k = 9 fc Sj, i.e., <9 fc Sj is a non-degenerate matrix. 

(iii) Vt a b = -d b E' a = 5° . Hence, E' a = -s a + D a (J x ). 

(iv) fij, = —d b E'\ i.e., S /l is affine in s a . 

In view of items (i) - (iv), the Liouville form E (|53|) reads 



E = x a dJ a + [D l (J x ) + BiiJjs^dJi + Si(Jj) V 

where we put 



a; a = -S ,a = s a -D a (J x ). (54) 



Since the matrix d Sj is non-degenerate, one can introduce new coordinates Ii = Sj(J. 
Ia — Ja- Then we have 

E = -x a dl a + [D^Ja) + B*(I x )s a \dIi + W. 

Finally, put 

^ = ^-[D'\h) + B'^h)s a } (55) 

in order to obtain the desired action-angle coordinates flTT|). These are bundle coordinates 
on U — > V where the coordinate shifts fl54|) - (^) correspond to a choice of another 



8 Appendix B 

In order to prove Proposition |3], we first show that functions i^F^ make up a CIS on the 
symplectic leaf (VqQ, Qq) and iVo is its invariant manifold without critical points (i.e., where 
first integrals fail to be independent). Clearly, the functions i^Fk are in involution, and A^" 
is their connected invariant manifold. Let us show that the set of critical points of {i^Fk} is 
nowhere dense in VqQ and iVo has none of these points. Let VqQ be equipped with some 
coordinates (<7 fc ,p;J. Then the trivial bundle £ flTj] ) is provided with the bundle coordinates 
(t, q k ,Pk) which play a role of the initial data coordinates on the momentum phase space V*Q. 
Written with respect to these coordinates, the first integrals F^ become time-independent. 
It follows that 

dF k (y) = df F k (Z(y)) (56) 

for any point y G V*Q. In particular, if y G VqQ is a critical point of {i^Fk}, then the 
trajectory ^ _1 (yo) is a critical set for the first integrals {Fk}. The desired statement at once 
follows from this result. 

Since N is compact and regular, there is an open neighbourhood of iV in V Q *Q isomorphic 
to V x N where V C M. m is a domain, and {v} x iVo, v G V, are also invariant manifolds in 
V *Q 0. Then 

W" = C\V xN ) = VxN (57) 

is an open neighbourhood in V*Q of the invariant manifold N foliated by invariant manifolds 
^ 1 ({ v } x ^o)> v G V, of the time-dependent CIS on V*Q. By virtue of the equality (|56|). 
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the first integrals {F k } have no critical points in W" . For any real number r 6 (— e, e), let 
us consider a section 

K : K*Q -> T*Q, p o fr r = g^) + r, 

of the affine bundle £ ©. Then the images /^(W 7 ") of W" ( |57|) make up an open neigh- 
bourhood U of h(N) in T*Q. Because ((U) = W", the pull-backs (*F k of first integrals 
Ffc are free from critical points in U, and so is the function 7i* dl2|) . Since the coordinate 
r = po — h provides a trivialization of the affine bundle (, the open neighbourhood U of h(N) 
is diffeomorphic to the product 

(-e,e) x ^ (-e,e) x V x h(iV) 

which is a trivialization of the fibration 

7i* x (xCF k ) : [/ -> (-e,e) x V. 

It remains to prove that the Hamiltonian vector fields of 7i* and C*F k on U are complete. 
It is readily observed that the Hamiltonian vector field (|13|) of 7i* is tangent to the 
manifolds h r (W" ), and is the image 77- = Th r o ^ H o ^ of the vector field 7^ fllOD. The latter 
is complete on W", and so is 7^ on [/. Similarly, the Hamiltonian vector field 

lk = -d t F k d° + d l FA - diF^ 

of the function (*F k on T*Q with respect to the Poisson bracket {, }t (§) is tangent to the 
manifolds h r (W"), and is the image 7^ = Th r o d k o ( of the Hamiltonian vector field d k of 
the first integral F k on W" with respect to the Poisson bracket {, }y The vector fields d k 
on W are vertical relative to the fibration W" — > R, and are tangent to compact manifolds. 
Therefore, they are complete, and so are the vector fields 7^ on {7. Thus, {/ is the desired 
open neighbourhood of the invariant manifold h(N). 
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